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A necessary and sufficient condition is given for the modulation- 
transfer- function of certain multimode optical fiber guides to be zero 
free in the closed right-half of the complex plane, and to be structurally 
stable with respect to that property. The condition is of interest, for 
example, in connection with the possibility of determining the phase 
of a modulation-transfer-function from its amplitude. 

I. INTRODUCTION AND PRELIMINARIES 

Reference 1 considers the range of validity of a Hilbert-transform 
approach in which the measured magnitude of the modulation-trans- 
fer-function of an optical fiber guide is used to compute the guide's im- 
pulse response.* It is argued there that a key "minimum-phase as- 
sumption" can fail to be satisfied in important cases, and a few closely 
related experimental and analytical results are presented. 

The purpose of this note is to report on a result along the same lines 
as a proposition given in Ref. 1 to the effect that, for a fiber guide that 
can propagate a finite number of discrete modes without mode mixing, 
the modulation-transfer-function (more precisely, the Laplace transform 
version of the modulation-transfer-function) is zero-free in the closed 
right half of the complex plane, and that property is structurally stable 
in a certain sense, if and only if a certain condition is met. The theorem 
described in Section II is concerned with a more realistic and far more 
interesting case in which mode mixing is not ruled out. In particular, the 
result provides further detailed support for the conclusion reached in 



* By "the guide's impulse response" is meant the output power of the guide excited by 
a unit impulse of optical power. The modulation-transfer-function G(w) is the envelope 
response of the fiber guide to an incoherent optical signal sinusoidally modulated at angular 
frequency w. To the extent that certain approximations hold, (Ref. 2), the impulse response 
is the Fourier transform of the modulation-transfer-function. The reason for considering 
the Hilbert transform approach is that it is often desirable to determine the impulse re- 
sponse of fiber guides by methods other than direct time-domain measurement, and, while 
|&(ui)| can be measured easily, it is at the present time difficult to accurately measure the 
phase of G(w). (See Refs. 3 and 4.) 
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Ref. 1. (i.e., that "nonminimum-phase behavior" is likely to arise, and 
can arise, in important actual cases). As in Ref. 1, it is the introduction 
and use of the concept of structural stability of a mathematical property 
which, for the case considered, enables a result to be obtained that is 
complete and easy to interpret.* 

We consider, as in Ref. 1, a class of optical fiber guides for which ap- 
propriate approximations can be made so that the modulation-trans- 
fer-function G (co) of a guide can be written in the form 



G(«) = £ 



e- iwT da(r), (1) 

T 

in which T denotes a closed, finite, real interval whose end points depend 
on the refractive indices of the core and cladding, and a(r) is a real-valued 
nondecreasing function.* It is assumed throughout that a(r) satisfies 
the normalization condition 



j; 



da{r) = 1. 

T 

As mentioned previously, in Ref. 1 attention is focused on the class of 
fiber guides that can propagate n discrete modes without mode mixing. 1 
In that case, G(co) can be written as 

n 

£ dje- iuT J, (2) 

in which each dj is a positive number that represents the initial excitation 
of the yth mode, and r\ < T2 < . . . < t„. 

In this note we suppose that the right side of (1) can be expressed 
as 

£ d,e-' w >+ C e' iuT b(r)dT, (3) 

;-l JT 

in which: hp is a positive integer (the motivation for using the subscript 
F will become clear shortly), d } > for 1 < j < kp, r\ < ti < . . . < r^ 
and 6(t) is a bounded piecewise-continuous § nonnegative function which 
takes into account mode mixing. It is also assumed that r\ has the fol- 
lowing property: 6(t) = for r e T with t < ti (which, of course, allows 
the possibility, but does not require, that t\ is equal to the lower endpoint 



* The general problem of determining when the Hilbert-transform approach (i.e., when 
the so-called Kramers-Kronig transformation) is valid is of interest in many fields (see, 
for example, Ref. 5). 

* Thus, roughly speaking, da(r) in (1) can be replaced with f(T)dr in which the function 
/(t) is nonnegative and may contain impulses corresponding to discrete modes. See Refs. 
2 and 3 for the relevant background material. We are assuming that material dispersion 
can be neglected. 

^Typically, n> 100. 

" The piecewise-continuity assumption appears to be adequate for applications. For 
basically a somewhat more general version of the theorem stated in Section II, see Section 
2.2. 
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of T). Since (3) is a specialization of the right side of (1), we also have 
Tj c T for each j , as well as 

Zdj+ f 6(T)dr=l. 
y-i Jt 

In physical terms, b is the relative-power density function associated 
with the nondiscrete modes, and the idealized impulse response of the 
guide is the inverse Fourier transform of (3). The observable impulse 
response of the guide (i.e., the impulse response of the guide-detector 
combination) is a somewhat smoothed version of the idealized response, 
with the smoothing provided by the detector (see Ref. 3). 

The important assumption that di > and that 6(t) = for t e T with 
t <t\ means that a discrete mode corresponding to the smallest modal 
delay is propagated. This assumption appears to be reasonable for at 
least some interesting classes of guides. For example, if a guide with a 
step-index profile is short enough to neglect mode conversion phe- 
nomena, then it is not unreasonable to assume that G(a>) has the form 
given in (2) with t\ the modal delay corresponding to the fundamental 
mode. In a real fiber, geometrical perturbations couple energy among 
the modes so that the distribution of modal delays changes continuously 
from a discrete set to a continuum as the fiber length L increases. Ex- 
perimental evidence indicates that the assumption is reasonable at least 
if the guide is not too long.* (For a particular fiber, there is a charac- 
teristic coupling length L c such that for L > L C) it is difficult in the time 
domain to isolate discrete modes with appreciable energy.) 

II. THE RESULT 

In this section, z denotes a complex variable and F(z) is defined by 
k 



F{z) = £ dje- ZT J+ f e- ZT b{T)d- 



(4) 



for each z. (Of course, if G(w) denotes (3), then G(u) = F(ia>), and F{z) 
is simply the Laplace transform of the generalized function whose 
Fourier transform is G(u)).) 

In order to state our result, consider an arbitrary function H(z) of the 
same type as F(z). More explicitly, let H(z) be given for all z by 

H(z)= £ 5;e-*<>+ ( e-"P{r)dT, (5) 

>=1 JT 

in which ku is a positive integer (not necessarily equal to kp), and the 
hj, the tj, and &(t) satisfy the restrictions imposed on the corresponding 
terms in (4). Let S denote the set of all such functions H{z). 



* The writer is indebted to his colleague I. P. Kaminow for a helpful discussion con- 
cerning the significance of the assumption described above. 
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For the purpose of defining a "distance" between F and an arbitrary 
element H of S, let Jf, Jh, and J, respectively, denote the sets of num- 
bers {1,2,. . . ,k F \, (1,2,. . . ,&//}, and {1,2,. . . ,min(^,/z//)}, and, with y a 
real variable, let q(y) denote any continuous nondecreasing function of 
y such that q(0) = 0. 

Let the "distance" p(F,H) between F and any H in S be defined 
by* 



(6) 



P(F,H)= E \dj-6j\+ E dj+ E *y 

;'£</ jt (Jf-J) j c (Jh-J) 

+ q (max \tj - tj\) + max I C° [6(r) - P(r)]d; 

Each term on the right side of (6) has a direct interpretation. In par- 
ticular, 

E dj+ E 6j, 

j c (Jf-J) j c (Jh-J) 

in which at most one sum is nonzero, reflects the extent to which terms 
in one of the two finite sums in (4) and (5) do not have counterparts in 
the other. Also, 



£ 



[b(r)-0(r))dr 



is an integral of the difference of two power-density functions, and, 
roughly speaking, if 



max 

u,u c T 



| £ [b(r) - 0(T)l<fc 



is sufficiently small, then, for practical purposes, the functions b and (3 
are indistinguishable in the sense that the observable impulse response 
of the guide is essentially unchanged if b is replaced with /?. (The portion 
of the idealized impulse response that does not contain impulses is g 
defined by g(t) = b(t) for t e T andg(t) = for t4 T. If, for example, the 
smoothing introduced by the detector is modeled by a filter with impulse 
response r given by r(t) = p~ l for t e [0 ,p] and r(t) = otherwise, in 
which p is a small positive constant, then the observable version of g(t) 
isp -1 ft- P g( T )dT for each t. Similarly, if instead r(t) = 0for£ <0,r(0) 
is finite, and the derivative of r is absolutely integrable on [0,°°), then 
an integration by parts shows that the observable version of g is essen- 
tially unchanged when b is replaced with a sufficiently nearby in the 
sense indicated above.) 

Our result, the theorem given below, provides an answer to the fol- 
lowing question: When is it true that F(z) ^ for Re(z) > and that 



* We adopt the convention that a sum over the empty set is zero. 
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property of F is structurally stable in the sense that there is a positive 
constant t (which can be thought of as a "tolerance") such that H(z) ^ 
for Re(z) > for every H in S such that p(F,H) < e. 
Theorem: We have F(z) ^ for Re(z) > with that property of F 
structurally stable, if and only if 



di> £ dt+ f b(r)dT. 



Jf 



Note that the theorem* does not rule out the possibility that the 
condition given in the theorem is violated and F(z) is zero-free in the 
closed right-half plane. (In fact, an example given in Ref. 1 shows that 
the possibility can occur. Essentially the same example can be used to 
show also that if the condition is violated, then it need not be the case 
that all functions in S "sufficiently close" to but different from F have 
a zero in Re(z) > 0.) On the other hand, at present it appears that there 
are complex, and for practical purposes impossible-to-specify, additional 
relationships among the tj, the dj, and b that, in particular take into 
account geometrical perturbations along the length of a real guide. The 
theorem shows that, when additional information is unavailable, it is 
not possible to prove that F(z) is zero-free in the closed right-half plane 
whenever 



di< £ dj + C b{r)dr 

i c Jf *sT 



j c Jf 



(which, in view of the normalization condition concerning a, is equivalent 
to the statement that the discrete mode with the smallest delay has 
relative power less than V 2 ) and, in the sense indicated above, the tj, the 
dj, and b are known only to within some tolerance, no matter how small 
the tolerance is. 

A proof of the theorem is given in the next section. 



* As indicated earlier, one application of the theorem is that it provides further detailed 
support for the material reported on in Ref. 1. For the benefit of the reader who has not 
read Ref. 1, we mention that a much simplified version of essentially the proof given in 
Section 2.1 can be used to show that if k F > 2, if 6(t) and /3(r) in (4) and (5), respectively, 
are each replaced by the zero function, if S is further restricted so that kn = kp and bj = 
dj for ;' = 1,2,. . . , kp for all H cS, and if p(F,H) is instead q(maxj c j\tj - tj \) [i.e., just 
the fourth term in the sum on the right side of (6)], then: F(z) ^ for Re(z) > and there 
is an i > such that H(z) ^ for Re(z) > for every H in the corresponding S with p(F,H) 
< e, if and only if d i > £y£ 2 dj. This result is basically a slight generalization of the com- 
parable proposition in Ref. 1. 
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2. 1 Proof of the theorem 

For the reader's convenience, we first repeat some of the material 
described above. We have T = [Ti,T 2 ] in which T\ and T 2 are real 
numbers such that T\ < T 2 , and S denotes the set of all functions of the 
complex variable z of the form (5) where kn is a positive integer, 8j > 
for 1 < ;' < k H , Ti < t x < t 2 < . . . < t kfi < T 2 , 0(r) is a nonnegative 
bounded piecewise-continuous* function defined on T such that (3{t) 
= for t e T with t < ri, and 

Z8j+ f/3(r)dr=l. 

7 = 1 JT 

The "distance" p(F,H) between any H cS and the particular element 
F of S given by (4), is defined by (6). 
Proof of the "If" Part: Suppose that 



d\ > £ dj + f b(T)dr. 
jcj F Jt 



if*i 



With 



j c Jf JT 



let c satisfy < e < (l/4)r. For each H e S with p(F,H) < e, we see that 
|di-5i| <€, 

E \dj-6j\ <e, E <*, <e, E 6; < e, 
;'c</ ;'c (Jf-J) j e (Jh-J) 

;V1 



as well as 



and therefore 



|J^(r)<iT- § T V{r)d 
h > E «y + f PWt. 



/f*i 



Thus, for Re(z) > and H c S with p(F,H) < e, we have 

\e zt *H{z)\ = 5i+ E 8je-*<*J-M + f e- 2 ^- f i)/3(r)dr 
j c jh Jt 



jvi 



* It will become evident that the theorem holds also if "piecewise continuous" is replaced 
with either "Riemann integrable" or "Lebesgue integrable." In this connection, see Sec- 
tion 2.2. 
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> 5i - E 5,e- z <0-'i) + f e-«< T -'»>j8(T)dT 

> «i - E «y - f j8(r)dT. (7) 

; c Jh J t 

Since the right side of (7) is positive, it is clear that H(z) 9± 0, which 

completes the proof of the "if" part. 

Proof of the "Only If" Part: Suppose now that 



di < E dj + f b{r)dr, 

jcj F Jt 



and let t > be given. 

Let k G = max(/2F,2), let 77 = min ((1/6) e, (1/2) di), and let 5 X = d 1 - 
77. If fc^ > 1, let 8 2 = d 2 + v and 8j = dj for ;' e \j:j e Jp'J 5^ 1,2}, and if fe^ 
= 1, let 82 = v and t 2 = T 2 . Then the function G defined (for all 2) by 

Kg r 

G{z)= E fy? -2T '"+ I c" 2T 6(T)dT 

belongs to S, and we have (if we set H = G): 

E K-5,|+ E dj+ E 5,-<^e (8) 

j eJ j c {Jf-J) i c (Jg-J) b 



and the strict inequality 

«1 < E «/+ f 6(r)dT. 



(9) 



;>i 



Let 5 T denote min{(r/+i - r,): 1 <;',; + 1 < fcc}, and let A = sup* c r 
6(t). With B = {t e T:b{r) > 0}, letsi and s 2 denote inf B and sup B, re- 
spectively, when B has nonzero measure, and T\ and T 2 , respectively, 
otherwise. 

Choose any 8 ( > such that g(5 e ) < (%) e, and let 5 be any positive 
number such that 

8 < min (\ 8 t , ] 8 T , \ e, e(18&)-\ \ (s 2 -sj). (10) 

\2 4 3 4 / 

Let u> = tr8~ l , and let K5 denote the set of numbers of the form t\ + k8, 
with k an odd positive integer. Clearly, exp[— iu)(t — ti)] = —1 for t e K&. 
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Choose t\ = t\, and (using kg ^ 2 and 2 5 < (%) 5 r ) for each; = 2,3,. . . ,&g 
choose a tj in K5 H [n.^VI such that | tj — tj\ < 25. Since 25 < (%) 5 T and 
25 < 5 e , we have fci < £2 < . . . < tk G and <?(max ; - e j\tj — tj\) < (%)e 

We see that the "distance" between F and the element E of S given 
by 

E(z) = E «/«"*'■' + f e- ZT 6(r)dT 

is at most (%) «. It therefore suffices to show that there is an H in S de- 
fined by 

H(z) = £ fye-**'"+ f e- 2T /3(r)dr 



with 



max I C [b(r) - P{t))<L 

u,v c T\ Ju 



1 
<-€ 

3 



such that H(z) = for some 2 with Re (2) > 0. 

Let L = K5 (si,S2)- Since 5 < Vt (S2 — si), 1 " L contains at least two 
points. Let the points in L be Pi,P2,- • • ,Pn, ordered so that Pi < P2 < 
. . . < p n . Let a be a positive number such that a <8,pi — s\> <r, and S2 
— p n > a. With I(u,u) denoting /{J b(r)dT for u,u e [si,S2], let /Mr) be 
defined for r e T by p\,( T ) = f(t - pi)/(s 1 ,p 1 + 5) + f(t - p 2 )I(p 2 - 5, p 2 
+ 5) + . . . + f(t - p n )I{p n - 5,s 2 ), in which /(t) = (2c-)- 1 for |t| < crand 
/(0 =0 otherwise. Since 



X 



/3„(r)dT = J(si,pi + 5)4- /(p 2 - 5,p 2 + 5) + . . . 

T 



we see that the function H a defined by 
ka 



+ I(p n -8,s 2 )= J T Hr)dr, (11) 
id by 
HM = E 8je-*tJ+ C PAr)e-"dT 

j=\ JT 

belongs to S. 

Using I(si,pi + 5) < 35A,* /(p, - 5,p y + 5) < 25A for; - 2, . . . ,(n - 
l),/(p n -5,s 2 )<35A, 

r f 6(r)dT= C PMdr 

for t = si,pi + 5,p2 + 5,. . . ,p n -\ + 5,s 2 , and the fact that b(r) and $ a (r) 
are nonnegative, we have 

I Cb{ T )dT- C P a {r)dr\ < 35A 



t See (10). 

* Notice that p\ — s\ < 25. 
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for all t cT.lt follows at once that for u,v c T, 

I JJ [b(r) - (3 a (T)]dr\ <65A 

3 

Therefore p(F,H a ) < e, uniformly for a as described. 
Let P(z) be defined for all z by 

P(z) = £ «,-e-»'J + e-*Pi/(si,pi + 5) 

7=1 

+ e" 2 P2/(p 2 - 5,p 2 + 5) + . . . + e-*P»I{p n ~ o,s 2 ). 

Let a be a real variable. Since tj c K& for; = 2,3,. . . ,kc, andpy c K b for 
j = 1,2,. . . ,n, and (9) and (11) hold, P(a + iu>) exp[(a + i«)tj* is real 
and negative when a = 0. On the other hand, P(a + ico) exp[(a + iu)ti] 
is positive for all sufficiently large a. Thus, P(zi) = for some z\ with 
Re(2,)>0. 

The function P(z) is analytic in z throughout the complex plane. Since 
it is not identically zero and is analytic at z — Z\, its zero at z = z\ is iso- 
lated. Therefore, there exists a constant r > such that r < Re(zi) 
and, with T = \z:\z - zi\ = r\, P(z) ^ for z c T. It follows that min 
{|P(z)|:z e r} is positive. 

Using the fact that 

• t+a 



(2a)- 1 f ° e~ ZT dT = e- zt w((jz), 

Jt-a 



in which w(az) = (2az)~ l (e za - e za ), we have 

HAz) = P(z) + [WW - l]y(z), (12) 

where 
y(z) = c-*^/(«uPi + o) + e-*P*I(p 2 ~ b,p 2 + o) 

+ . . . + e-*P»I{p n -8,s 2 ). 

The function \y(z) |is bounded on T, and max \\w(az) — l\:z c T\ (and 
hence max }|[u;((tz) - l]y(z)|:z c T\) can be made arbitrarily small by 
choosing a to be sufficiently small. By Rouche's theorem, for sufficiently 
small o, H a {z) has a zero inside T (and hence in Re(z) > 0). This com- 
pletes the proof. 



♦ See (10). 

1 Recall that w = rf" 1 . 
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2.2 Comment 

It is apparent that we have proved a somewhat stronger result than 
the theorem stated. Suppose that the definition of S is changed to the 
extent that /3(r) need not be piecewise continuous, but merely Lebesgue 
integrable. Then it it is clear that the "if part of the theorem remains 
true. More importantly, the proof shows that if 



di< Z d J + f b(r)dT, 



J cJf 

;Vl 

then, for any e > 0, there is an H in S of the form (5) with the following 
properties: p(F,H) < e, H(z) has a zero in the open right-half of the 
complex plane, t\ = r\,hn = kp t bj = dj for; = 1,2,. . . ,kp, (3{t) is piecewise 
constant, and the smallest closed real interval containing the support 
of 6(t) (which might possibly be the "empty interval") also contains the 
support of 0(t). 
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